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Abstract 

We study optimal algorithms in adaptive sampling recovery of smooth functions defined 
on the unit d-cubc I'' := [0, 1]'^. The recovery error is measured in the quasi-norm || • ||q of 
Lq := Lq{I'^). For B a subset in Lq, we define a sampling recovery algorithm with the free 
choice of sample points and recovering functions from B as follows. For each / from the quasi- 
normed Besov space B^g, we choose n sample points. This choice defines n sampled values. 
Based on these sample points and sampled values, we choose a function from B for recovering 
/. The choice of n sample points and a recovering function from B for each / G B^g defines 
a n-sampling algorithm by functions in B. We suggest a new approach to investigate the 
optimal adaptive sampling recovery by in the sense of continuous non-linear n- widths which 
is related to ri-term approximation. If $ = {y^k}keK is a family of elements in Lg, let !]„($) 
be the non-linear set of linear combinations of n free terms from $, that is E„($) := {ip = 
Sj=i '^j'Pkj '■ kj ^ K}. Denote by Q the set of all families $ in Lg such that the intersection 
of $ with any finite dimensional subspace in Lg is a finite set, and by C{Bpg, Lg) the set of all 
continuous mappings from B^g into Lg. We define the quantity 

i/„(5"e,L„) := inf inf sup \\f-S^{f)\\g. 

Let < p,q^9 < oo and a > d/p. Then we prove the asymptotic order 

We also obtained the asymptotic order of quantities of optimal recovery by in terms of best 
n-term approximation as well of other non-linear n-widths. 

Keywords Adaptive sampling recovery • n-sampling algorithm • B-spline quasi-interpolant 
representation • B-spline • Besov space 

Mathematics Subject Classifications (2000) 41A46 • 41A05 • 41A25 ■ 42C40 
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1 Introduction 



The purpose of the present paper is to investigate optimal algorithms in adaptive sampling recovery 
of functions defined on the unit d-cube I'^ := [0, l]'^. Let Lq := Lq{I'^), < q < oo, denote the 
quasi- normed space of functions on I'^ with the usual gth integral quasi- norm || • \\q for < g < oo, 
and the normed space C(I'^) of continuous functions on with the max-norm || • ||oo for q = oo. 
For < p,6,q < oo and a > 0, let B^q be the quasi-normed Besov space with smoothness a, 
equipped with the quasi- norm || • \\b" (see Section [2] for the definition). We consider problems 
of adaptive sampling recovery of functions from B^g. The recovery error will be measured in the 
quasi- norm || • ||g. 

We first recall some well-known non-adaptive sampling recovery algorithms. Let X be a quasi- 
normed space of functions defined on W^, such that the linear functionals / i— ?> f{x) are continuous 
for any x G I'^. We assume that X C Lq and the embedding Id : X ^ Lg is continuous, where 
Id(/) := /. Suppose that / is a function in X and ^„ = {a^'^}fc=i are n points in I"^. We want to 
approximately recover / from the sampled values f{x^), f{x'^), fix"")- A classical linear sampling 
algorithm of recovery is 

n 

Ln{f) = Ln{^n,inJ) ■■= ^f{x'')^k, (LI) 

k=l 

where = {fk}k=i are given n functions in Lq. A more general sampling algorithm of recovery 
can be defined as 

Rn{f) = Rn{Hn,CnJ) := HM (x^) , f {x'')) , (L2) 

where Hn is a given mapping from to Lq. Such a sampling algorithm is, in general, non- linear. 
Denote by SX the unit ball in the quasi-normed space X. To study optimal sampling algorithms 
of recovery for f & X from n their values by algorithms of the form (jl.2p , one can use the quantity 

9n{X,Lq) := inf sup \\f - Rn{Hn,^n, f)\\q, 

where the infimum is taken over all sequences = {x'^}'^^^ and all mappings Hn from M" into Lq. 

We use the notations: x+ := max(0,x) for x e M; An{f) < Bn{f) if An{f) < CBn{f) with C 
an absolute constant not depending on n and/or f € W, and An{f) x Bn{f) if An{f) <C Bn{f) 
and An{f) <C An{f). It is known the following result (see [13], [23], [26], [27], [31] and references 
there). If < p,9,q < oo and a > d/p, then there is a linear sampling recovery method L* of the 
form (jl.ip such that 

gn{B^,e,Lq) x sup ||/-L;(/)||, x n--/'^+(i/f-i/9)+. (1.3) 

This result says that the linear sampling algorithm L* is asymptotically optimal in the sense that 
any sampling algorithm i?„ of the form (|1.2p does not give the rate of convergence better than L* . 

Sampling algorithms of recovery of the form (jl.2p which may be linear or non-linear are non- 
adaptive, i.e., the points ^„ = {x'^}^^^ at which the values f{x^),...,f{x"') are sampled, and the 
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method for construction of recovering functions are the same for ah functions / G X. Let us 
introduce a setting of adaptive samphng recovery. 

If i? is a subset in Lg, we define a samphng algorithms of recovery with the free choice of 
sample points and recovering functions from B as follows. For each / € X we choose a set of n 
sample points. This choice defines a collection of n sampled values. Based on the information of 
these sampled values, we choose a function from B for recovering /. The choice of n sample points 
and a recovering function from B for each / G X defines a sampling algorithms of recovery by 
functions in B. More precisely, a formal definition of is given as follows. Denote by Z" the set 
of subsets ^ in of cardinality at most n, V" the set of subsets r/ in M x 1°' of cardinality at most n. 
Let Tn be a mapping from X into I". Then T„ generates an n-sampling operator /„ from X into 
V" which is defined as follows. If T„(/) = then /:„(/) = {(/(x^), x^), (/(x"), x'^)}. 

Let a mapping from V" into B. Then the pair {In,Pn) generates the mapping from X 
into B, by the formula 

S^{f):=P^{lM)), (1-4) 

which defines a n-sampling algorithm with the free choice of n sample points and approximant 
from B. We call the mapping a recovering operator. 

Clearly, a linear sampling algorithm L„(<^„,^„,-) defined in (|l.ip is a particular case of 5^. 
We are interested in adaptive n-sampling algorithms of special form which are an extension of 
Ln{^n,^n, ") to an n-sampling algorithm with the free choice of n sample points and n functions 

= {v^k}k=i for each / € X. To this end we let $ = {'fk}keK be a family of elements in Lg, and 
consider the non-linear set S„(<I>) of linear combinations of n free terms from $, that is 

n 

:= { 99 = ^ ajifkj : kj £ K}. 
i=i 

Then for B = an n-sampling algorithm is of the following form 

keQiv) 

where (pk € ^, rj = In{f), Qiv) K with |Q(?7)| < n and are functions on V". 

We want to choose an n-sampling algorithm so that the error of this recovery ||/ — S^{f)\\q 
is as smaller as possible. Clearly, such an efficient choice should be adaptive to /. To investigate 
the optimality of (non-continuous) adaptive recovery of functions / from the quasi-normed space X 
by n-sampling algorithms of the form (|1.5p . we introduce the quantity as follows: 

Sn{X,^,Lg) := inf sup ||/ - 5^(/)||,. (1.6) 

The definition ()1.6p corrects a definition of s„(X, Lg) which has been introduced and denoted 
by i'n{SX,^)g and s„(5X, $)g in [l^ and [T8], respectively. The quantity s„(X, <I>,Lg) is directly 
related to non-linear n-term approximation. We refer the reader to [32j for surveys on various 
aspects in the last direction. 
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The quantity ^,Lq) depends on the family <I> and therefore, is not absolute in the sense 

of n- widths or optimal algorithms. We suggest an approach to investigate the optimal adaptive 
sampling recovery by in the sense of continuous non-linear re-widths which is related to re-term 
approximation too. Namely, we consider the optimality in the restriction with only n-sampling 
algorithms of recovery of the form (jl.Sp and with a continuity assumption on them. Continuity 
assumptions on approximation and recovery algorithms have their origin in the very old Alexandroff 
re-width which characterizes best continuous approximation algorithm by re-dimensional topological 
complexes [Ij (see also [33] for details). Later on, (continuous) manifold n- width was introduced 
by in ^24j, and investigated in [12], [9], [20], [H], [I5],[l6]. Several continuous re-widths based 
on continuous algorithms of re-term approximation, were introduced and studied in [H], [T5].|16|. 
The continuity assumption is quite natural: the closer objects are the closer their reconstructions 
should be. A first look seems that a continuity restriction may decrease the choice of approximants. 
However, in most cases it does not weaken the rate of the corresponding approximation. Continuous 
and non-continuous algorithms of nonlinear approximation give the same asymptotic order. This 
motivate us to consider continuous n-sampling algorithms of recovery . Since we assume that 
functions to be recovered are living in the quasi-normed space X and the recovery error is measured 
in the quasi-normed space Lg, the requirement that G C{X, Lq) is quite reasonable. (Here and 
in what follows, C{X,Y) denotes the set of all continuous mappings from X into Y for the quasi- 
metric spaces X,Y). This leads to the following definition. 

Denote by G the set of all families $ in Lq such that the intersection of $ with any finite 
dimensional subspace in Lq is a finite set. We define the quantity 

Un{X,Lq) := inf inf sup \\f-S^{f)\\q. 

The restriction <I> G ^ in the definition of Vn{X,Lq) is minimal and natural for all well-known 
approximation systems. 

Another way to study optimal adaptive (non-continuous) re-sampling algorithms of recovery 
in the sense of nonlinear re- widths has been proposed in |17j . |18j . In this approach, B is required 
to have a finite capacity which is measured by their cardinality or pseudo-dimension. Given a 
family B of subsets in Lg, we consider optimal sampling recoveries by B from B in terms of the 
quantity 

Rn{W,B)q:= inf inf sup ||/ - S^(/)||,. (1.7) 

-See SB f(zw 

We assume a restriction on the sets B ^ B, requiring that they should have, in some sense, a 
finite capacity. In the present paper, the capacity of B is measured by its cardinality or pseudo- 
dimension. This reasonable restriction would provide nontrivial lower bounds of asymptotic order 
of Rn{W,B)q for weh known function classes W . Denote Rn{W,B)q by en{W)q if B in (jl.7p is 
the family of all subsets B in Lq such that \B\ < 2", where \B\ denotes the cardinality of B, and 
by rn{W)q if B in ()1.7p is the family of all subsets B in Lq of pseudo-dimension at most re. The 
definition (jl.7p corrects definitions of en{W)q and rn(W)q introduced in [18] . 

The quantity e„(VF)q is related to the entropy re-width (entropy number) en{W)q which is the 
functional inverse of the classical e-entropy introduced by Kolmogorov and Tikhomirov [22] . The 
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quantity rn{W)g is related to the non-linear n- width pn{W)q introduced recently by Ratsaby and 
Maiorov [29]. (See the definition of en{W)q and pn{W)q in Section [5]). 

The pseudo-dimension of a set B of real- valued functions on a set 0, is defined as follows. 
For a real number t, let sgn(t) be 1 for t > and —1 otherwise. For x G M", let sgn(x) = 
(sgn(xi), sgn(x2), sgn(x„)). The pseudo-dimension of B is defined as the largest integer n such 
that there exist points a} , . . . in Q. and b G R" such that the cardinality of the set 

{ sgn(y) : y = {f{a') + h, f{a^) + &2, • • • , /(«") + &n), / G 5} 

is 2". If n is arbitrarily large, then the the pseudo-dimension of B is infinite. Denote the pseudo- 
dimension of B by dimp(i?). The notion of pseudo-dimension was introduced by Pollard [28] and 
later Haussler [21] as an extention of the VC-dimension suggested by Vapnik-Chervonekis for 
sets of indicator functions. The pseudo-dimension and VC-dimension measure the capacity of a 
set of functions and are related to its e-entropy (see also [2U], ^30j). If S is a n-dimensional linear 
manifold of real- valued functions on fi, then dimp(i?) = n (see |21j). 

We say that p, q, 0, a satisfy Condition (jl.Sp if 

< p, q,0 < oo, < a < oo, and there holds one of the following restrictions : 

(i) a > d/p; (1.8) 

(ii) a = d/p, 9 < min(l,p), p,q < oo. 

Let M be the set of B-splines which are the tensor product of integer translated dilations of the 
centered cardinal spline of order 2r (see the definition in Section [2|). 

The main results of the present paper are read as follows. 

Theorem 1.1 Let p,q,9,a satisfy the Condition (II. 8p and a < 2r. Then for the d-variable Besov 
space B^Q, there is the following asymptotic order 

s„(i?°„M,L,) X MB'^,e,Lq) x r„(5i?°,)g x e„(5i?°,), x n""/-^. (1.9) 

Comparing this asymptotic order with (jl.3p . we can see that for < p < q < oo, the asymptotic 
order of optimal adaptive sampling recovery in terms of the quantities Sn, Vn-, and r„, is better 
than the asymptotic order of any non-adaptive n-sampling algorithm of recovery of the form ()1.2p . 

To prove the upper bound for (jl.9p . we use a B-spline quasi-interpolant representation of 
functions in the Besov space B^q associated with some equivalent discrete quasi-norm [T7], [18j . 
On the basis of this representation we construct corresponding asymptotically optimal n-sampling 
algorithms of recovery which give the upper bound for (jl.9p . The lower bound of ()1.9p is established 
by the lower estimating of the smaller related n-widths and the quantity of n-term approximation. 

The paper is organized as follows. 

In Section [21 we give a definition of quasi-interpolant for functions on 1°', describe a B-spine 
quasi-interpolant representation for Besov spaces B^^ with a discrete quasi-norm in terms of the 
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coefficient functionals. The proof of the asymptotic order of Vn[B'^Q,Lq) in Theorem 11.11 is given 
in Sections [3] and [H More precisely, in Section [3l we construct asymptotically optimal adaptive 
n-sampling algorithms of recovery which give the upper bound for i'n{Bpg, Lg) (Theorem 13. ip . In 
Section [4] we prove the lower bound for L^) (Theorem 14. ip . In Section [H we prove the 

asymptotic order of SniB^g,M, Lg), rn{SB^g)g and eniSB^g)g in Theorem[LTl 



2 Preliminary background 

For a given natural number r, let M be the centered B-spline of even order 2r with support [— r, r] 
and knots at the integer points — r, ...,0, ...,r and define the B-spline 

Mk,six) := M(2^x-s), 

for a non-negative integer k and s € Z. To get the d-variable B-spline M^^g for a non-negative 
integer k and s € Z*^, we let 

d 

M{x) := JjM(xi), X = (xi,X2, 

i=l 

and 

Mk,s{x) := M{2''x-s). 
Denote by M the set of all ,, which do not vanish identically on I'^. 

Let A = {A(j)}jgpd(^) be a finite even sequence, i.e., A(— j) = A(j), where -P'^(/u) := {j G Z'^ : 
^ ^ = 1; 2, d}. We define the linear operator Q for functions / on M'^ by 

Q(/,x):= A(/,s)M(x-s), (2.1) 

where 

A(/,.):= Yl ^U)f('-j)- (2-2) 
The operator Q is bounded in C(IR'^) and 

IIW)llc(M'^)<l|A||||/|lc(R-) 

for each / e C(M'^), where 

iiAii= Yl 

j€P''{^i) 

Moreover, Q is local in the following sense. There is a positive number 5 > such that for any 
/ € C(M'^) and x € W^, Q{f, x) depends only on the value f{y) at a finite number of points y with 
lUi — Xi\ < 6, i = 1, 2, ...d. We will require Q to reproduce the space V2r-i of polynomials of order 
at most 2r — 1 in each variable Xi, that is, 

Q{p) = p, p^K^v 
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An operator Q of the form ()2.ip - (|2.2p reproducing "P^r-i' called a quasi-interpolant in C{M.'^). 

There are many ways to construct quasi-interpolants. A method of construction via Neumann 
series was suggested by Chui and Diamond [4j (see also [3l p. 100-109]). De Bore and Fix 
[S] introduced another quasi-interpolant based on the values of derivatives. The reader can see 
also the books [3], for surveys on quasi-interpolants. The most important cases of d-variate 
quasi-interpolants Q are those where the functional A is the tensor product of such d univariate 
functionals. Let us give some examples of univariate quasi-interpolants. The simplest example is 
a piecewise linear quasi-interpolant is defined for r = 1 by 

Q{f,x) := ^f{s)M{x-s), 

sez 

where M is the symmetric piecewise linear B-spline with support [—1, 1] and knots at the integer 
points —1,0,1. This quasi-interpolant is also called nodal and directly related to the classical 
Faber-Schauder basis p^. Another example is the cubic quasi-interpolant defined for r = 2 by 

Q{f,x) := ^l{-/(s-l) + 8/(5)-/(s + l)}M(x-s), 

where M is the symmetric cubic B-spline with support [—2, 2] and knots at the integer points 
-2,-1,0,1,2. 

Let Q, = [a,b]'^ be a d-cube in U.'^. Denote by Lp{Q) the quasi-normed space of functions on 
17 with the usual pth integral quasi-norm || • ||p^Q for < p < oo, and the normed space C{Q) of 
continuous functions on J7 with the max-norm || • \\oo,n for p = oo. If r be a number such that 
< r < min(p, 1), then for any sequence of functions {fk} there is the inequality 

\\T.f4ln - T.Wf^'Wl^- (2.3) 



We introduce Besov spaces of smooth functions and give necessary knowledge of them. The 
reader can read this and more details about Besov spaces in the books [2], [25], [10]. Let 

w«(/,i)p := sup ||A^/||p^id(,;,) 

\h\<t 

be the Ith modulus of smoothness of / where l'^{lh) := {x : x,x + Ih & I*^}, and the lih. difference 
A^/ is defined by 



i=o ^-^^ 



f{x + 3h). 



For < p, 6* < oo and < a < /, the Besov space B^q is the set of functions f G Lp for which the 
Besov quasi-semi-norm I/Ibq^ is finite. The Besov quasi-semi-norm |/|_B"g is given by 

\sup,>o t-"^z(/,i)p, = ^- 
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The Besov quasi-norm is defined by 

Bif) = \\f\\B^y-= + 



If {/fclfcLo ^ sequence whose component functions fk are in Lp, for < p, < oo and f3 > 
we use the hQ{Lp) "quasi-norms" 



with the usual change to a supremum when = oo. When {/fej^g ^ positive sequence, we 
replace ||/fc||p by \fk\ and denote the corresponding quasi-norm by ||{/fe}|L/3. 

For the Besov space B'^q, there is the following quasi-norm equivalence 
B{f) X := \\{i^i{f\2-%]\\,<. + ll/llp. 



In the present paper, we study the sampling recovery of functions from the Besov space B^q 
with some restriction on the smoothness a. Namely, we assume that a > d/p. This inequality 
provides the compact embedding of B'^g into C(I'^). In addition, we also consider the restriction 
a = d/p and 9 < min(l,p) which is a sufficient condition for the continuous embedding of B^q 
into C{I'^). In both these cases, B!^g can be considered as a subset in C(I'^). 

If Q of is a quasi-interpolant of the form (12.ip - (l2.2p . for h > and a function / on M'^, we 
define the operator Q{-;h) by 

Q(.f;h) = a-hoQ oai/hif), 
where ah{f,x) = f{x/h). By definition it is easy to see that 

Qif,x;h)= ^Aif,k;h)Mih-'x-k), 
k 

where 

A{f,k;h):= K3)f{h{k-3)). 
j€P''{^i) 

The operator Q{-;h) has the same properties as Q: it is a local bounded linear operator in M"^ 
and reproduces the polynomials from V2r-i- Moreover, it gives a good approximation of smooth 
functions [6, p. 63-65]. We will also call it a quasi-interpolant for C(M'^). 

The quasi-interpolant h) is not defined for a function / on I'^, and therefore, not appropriate 
for an approximate sampling recovery of / from its sampled values at points in I'^. An approach 
to construct a quasi-interpolant for a function on I'^ is to extend it by interpolation Lagrange 
polynomials. This approach has been proposed in [17] for the univariate case. Let us recall it. 

For a non-negative integer m, we put Xj = j2~^,j G Z. If / is a function on I, let Umif) and 
Vm be the (2r — l)th Lagrange polynomials interpolating / at the 2r left end points xq, xi, ...,X2r-i, 
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and 2r right end points 3;2'"-2r+ii 3;2™-2r+3) 2:2™ , of the interval I, respectively. The function 
fm is defined as an extension of / on M by the formula 



fmix) := < 



Um{f,x), X<0, 
fix), 0<X<1, 
yrn{f,x), X>1. 



Obviously, if / is continuous on I, then fm is a continuous function on M. Let Q be a quasi- 
interpolant of the form (j2.ip - (j2.2p in C(M). We introduce the operator Qm by putting 

Qrn{f,x) := Qifm,x;2-^), X € I, 

for a function / on I. By definition we have 

Qm{f,x) = ^ am,s{f)Mm,s{x), \/x el, 

where J(m) := {s G Z : —r < s < 2™ + r} is the set of s for which Mm,s do not vanish identically 
on I, and 

amAf)-= A(/m,5;2-")= ^ A(j)/m(2-"^(5-j)). 
The multivariate operator is defined for functions / on by 

sGJ{m) 

where J'^{m) := {s G Z*^ : — r < Sj < 2™ + r, i = 0, 1, ...,d} is the set of s for which Mm,s do not 
vanish identically on f^, and 

if) = am,si{{a 

m,S2 [■ ■ ■^m,s^ if))), (2.4) 

where the univariate functional Om.Si is applied to the univariate function / by considering / as a 
function of variable Xi with the other variables held fixed. Moreover, the number of the terms in 
Qmif) is of the size ^ 2'^"^. 

The operator Qm is a local bounded linear mapping in (7(1'^) and reproducing "P^r-i- P^^' 
ticular, 

||Q™(/)||c(i^)<C||A||||/||c(i.) (2.5) 
for each / G C(I"'), with a constant C not depending on m, and, 

Qmip*) = P, P(^K-1, (2-6) 

where p* is the restriction of p on W^. The multivariate operator Qm is called a quasi-interpolant 
in C(I^). 



From (j2.5p and (j2.6p we can see that 

ll/-Qm(/)lb(irf) ^0, m^oo. (2.7) 
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Put M(m) := {Mm,s G M : s G J"'(m)} and V(m) := spanM(m). U < p < oo, for all 
non-negative integers m and all functions 

g= ('sMm,s (2.8) 

from V(m), there is the norm equivalence 

\\g\\p X 2-'^-/P||{aJ||p,^, (2.9) 

where 



i/p 



with the corresponding change when p = oo. 

For non-negative integer A:, let the operator qk be defined by 

%(/) :=Qfc(/)-Qfc-i(/) with Q_i(/):=0. 
From (12. 6p and (12. 7p it is easy to see that a continuous function / has the decomposition 



k=0 

with the convergence in the norm of C{I'^). By using the B-spline refinement equation, one can 
represent the component functions Qkif) as 

(ik{f)= Yl (^kAfWk,s, (2.10) 

where c^^s are certain coefficient functionals of /, which are defined as follows. For the univariate 
case, we put 

CkAf) ■= «m(/) - ^ > 0, (2.11) 

4,,(/) := 2-2'-+! Y C^')ak-i,m{f), k>0, a^,,(/) := 0. 
{m,j)eC{k,s) ^ ^ 

and 

C{k,s) := {{m,j) ■.2m + j - r = s, me J{k-l), < j < 2r}, A; > 0, C(0,s) := {0}. 
For the multivariate case, we define Ck,s in the manner of the definition (j2.4p by 

Ck,s{f) '■= Ck,si{{ck,S2{-Ck,sa{f)))- (2-12) 

For functions / on I'^, we introduce the quasi- norms: 
B2if) ■■= \\{qk{f)}\\b^iL,y, 

i/e 

B,if) := ( j;(2("-'^/^)^-||{c,,,(/)}||,,,)'" 
The following theorem has been proven in [18] . 
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Theorem 2.1 Let < p,9 < oo and d/p < a < 2r. Then the hold the following assertions, 
(i) A function f € B^q can be represented by the mixed B-spline series 

oo oo 

/ = E ^kif) = E E ^kAf)Mk,s, (2.13) 

k=0 fc=OsGJ<*(fe) 

satisfying the convergence condition 

B2{f) - Bsif) « B{f), 

where the coefficient functionals Ck^sif) o,re explicitly constructed by formula (|2.11|) - (j2.12p as 
linear combinations of at most N function values of f for some N gN which is independent 
of k, s and f . 

(a) If in addition, a < min(2r, 2r — 1 + 1/p), then a continuous function f on I"^ belongs to the 
Besov space B^g if and only if f can be represented by the series (j2.13p . Moreover, the Besov 
quasi-norm B(f) is equivalent to one of the quasi-norms B2{f) and B^{f). 

3 Adaptive continuous sampling recovery 

In this section, we construct asymptotically optimal algorithms which give the upper bound of 
Vn{BpQ,Lq) in Theorem ll.il We need some auxiliary lemmas. 

Lemma 3.1 Let p,q,0,a satisfy Condition (jl.Sp and a < 2r. Then Qm € C{B'^Q,Lq) and for any 
f G BpQ, we have 

||Qm(/)||, < C||/||b.^, 

< C'2-^-'iWp~yi)+y'-\\f\\s^^^ (3.1) 
with some constants C,C' depending at most on d,r,p,q,9 and ||A||. 

Proof. We first prove (j3.ip . The case when the Condition (ii) holds has been proven in [T5]. Let 
us prove the case when the Condition (i) takes place. We put a' := a — d{l/p — ^/q)+ > 0. For 
an arbitrary / G B'^q, by the representation (I2.13P and (|2.3p we have 

\\f-Qm{f)\\l « E (3.2) 
k>m 

with any r < min(g, 1). From ()2.10p and ()2.8|) - ()2.9p we derive that 

||g,(/)||, « 2(i/P-V.)+fc||^^(/)||^ (3.3) 
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Therefore, if < mm{q, 1), then we get 




« 2-"'-i?2(/) « 2--'-||/||b"_^. 
Further, ii 6 > inm{q, 1), then from (I3.2p and (I3.3p it fohows that 

\\f-QM\\( « T.h>^(f)\\( 

k>m 

« Y.{2-%k{f)\ur'{2~'''''r\ 

k>m 

where q* = min(g, 1). Putting u := 9/q* and i^' := v jiy — 1), by Holder's inequahty obtain 

\\!-Qmmt « (Ei2"'iift(/)ii,r*^] (Ei2""''}'*'^') 
« {B2(/)r*{2-"'-r* « {2-"'-ri/iii;. 

Thus, the inequahty (IS.ip is completely proven. 
Next, by use of the inequality 

with any r < min(g,l), in a similar way we can prove that ||(5m(/)||ij ^ ll/ll-B°g ^"^^ therefore, 
the inclusion Qm G ^i^pe^ ^g)- '-' 

Lemma 3.2 For functions f on Y^, Qk defines a linear n-sampling algorithm of the form (jl.ip . 
More precisely, 

sG/'*(fc) 

where n := {2^ + 1)"^, j are explicitly constructed as linear combinations of at most {2fi + 2)'^ 
B-splines Mk,s, and I'^{k) := {s G : < Si < 2^=, i = 1, d}. 
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Proof. For univariate functions the coefficient functionals ak,s{f) can be rewritten as 
akAf) = E ^(^-i)A(2~'i) = E ^M(j)/(2-'i), 

\s-j\<n jeP{k,s) 

where Afc,,(j) := A(s - j) and P(A;, s) = := {j G {0, 2*^} : s - j G P(^)} for ^ < s < 2*^ - 

^k,sU) is a hnear combination of no more than max(2r, 2/x + 1) < 2/Lt + 2 coefficients X{j),j G -P(/x), 
for s < or s > 2*^ — and 



P,(/i)U{0,2r-l}, s</x, 
Ps(m)U{2'= -2r + l,2'=}, s>2^-/x. 



If J G P{k,s), we have |j — s| < max(2r, /i) < 2/x + 2. Therefore, P{k,s) C PsifJ-), and we can 
rewrite the coefficient functionals 0^,5 (/) in the form 

afe,s(/) = E >^kM)m~'j) 

j-seP(2n+2) 

with zero coefficients Xk,s{j) for j ^ P{k, s). Therefore, for any k G we have 



r 

s 



QkU) = E ^kAmi,s = E E A,,,(i)/(2-^i)M,^ 

seJ(A:) seJr{k)j-seP{2fj,+2) 

= E E Tfcj(^)M,^,, 

je/(fe) s-jeP(2(U+2) 

for certain coefficients jk,j{s)- Thus, the univariate qk{f) is of the form 
where 

V'feJ := E ^k,j{s)Mk,s, 

s-j6P(2//+2) 

arc a hncar combination of no more than the absolute number 2)U + 2 of B-splines M^^s^ and the 
size \I{k)\ is 2^. Hence, the multivariate qk{f) is of the form 

Qkif) = E /(2~'i)V'ikj, 

where 

d 

V'fcj := Yli'kJi 
1=1 

are a linear combination of no more than the absolute number (2/x + 2)*^ of B-splines 5, and 
the size \I'^{k)\ is 2'^''. □ 
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For < p < oo, denote by the space of all sequences x = {x^}^]^ of numbers, equipped 
with the quasi-norm 

with the change to the max norm when p = oo. Denote by the unit ball in i^. Let £ = {ek}^=i 
be the canonical basis in i. e., x = YlT=i ^k^k- 

For X = {xk}^^i € i"^, we let the set {kj}'jLi be ordered so that 

We define the algorithm P„ for the n-term approximation with regard to the basis £ in the 
space £^ (n < in) as follows. For x = {xk}^^i G we let the set {kj}J!^i be ordered so that 

\ > IXjjl > • • • \^js \ ^ ■ ■ ■ ^ I^^Jrnl- 

Then, for n < m we define 

n 

Pn{x) := "^{xk^ - |x„+i|sign2;fc^.)efc.. 
For a proof of the following lemma see |16) . 

Lemma 3.3 The operator Pn € C{P^, l^) for < p,q < oo. If < p < q < oo, then we have for 
any positive integer n < m 

sup \\x-Pn{x)\\ern < n^/'?-Vp. 



The following theorem gives the upper bound of ()1.9p in Theorem 11.1 

leorem 3.1 Let p,q,6,a satisj 
there is the following upper bound 



Theorem 3.1 Let p,q,6,a satisfy Condition (jl.Sp . Then for the d-variable Besov space B^q, 



UniB^e^L,) « n-^l". (3.4) 

If in addition, a < 2r, we can find an positive integer k* and a continuous n-sampling recovery 
algorithm G C(5^g, Lq) of the form (jLll) with A = S„(M(A:*)), such that 

sup ||/-5^(/)|l,«n-"/'^. (3.5) 



Proof. We will prove (13. 5p and therefore, (|3.4p . We first consider the case p > q. For any integer 
k* , by Lemmas 13.21 and 13.11 we have 

sup \\f-Qk*{f)h - 2-"*^*- (3-6) 
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The number of sampled values of / in Qk* (/) is n* := (2'^* + !)'^. For a given integer n (not smaller 
than 2'^), define k* by the condition 



Cn<n* = (2'=* + if < n, (3.7) 

with C an absolute constant. By Lemma [3?T] Qi,* € C{B'^ q, Lg). By the choice of k*, Qk*{f) = 
S^{f) is a linear n-sampling algorithm S^{f) is of the form p.ip with A = S„(M(A;*)) and 
M(A;*) € as a finite family. Therefore, by (j3.6p and (j3.7p we receive (jS.Sp for the case p > q. 

We next treat the case p < q. For arbitrary positive integer m, a function / G SB^q can be 
represented by a series 

m oo 

/ = E E 2^ 2^ Cfc,,(/)Mfc,, (3.8) 

converging in the norm of B^g or, equivalently, 

oo 

/ = Qm{f)+ "Ikif) (3.9) 

fc=m+l 

with the component functions 

Qk{f)= Yl ^kAf)Mk,s (3.10) 

SgJ<*(fc) 

from the subspace V(A;). Moreover, q'fc(/) satisfy the condition 

\\qkif)\\p >i 2-*/f||K,(/)}||p,fc « 2-"^ A: = m + l,m + 2,... (3.11) 

The representation (|3.8p - (j3.1ip follows from Theorem 12.11 for the case (i) in Condition (jl.Sp . and 
from Lemma |3. II for the case (ii) in Condition (jl.Sp . 



Put m{k) := \ J'^{k)\ = {2^ + 2r — I)"'. Let k,k* be non-negative integers with k < k* , and 
{n(fc)}^^^_i_-^ a sequence of non-negative integers with n{k) < m{k). We will construct a recovering 
function of the form 

k" n{k) 

G{f):= <^kAfWk,s + Y Y''k,s,{f)Mk,s,, (3.12) 

seJ(fc) k=k+i i=i 

with Sfc j G J'^{k), or equivalently, 

G(/) = Q-kif) + Y ^'^(fy (3.13) 

k=k+l 

The algorithms G/; are constructed as follows. For a / G SB^g, we take the sequence of 
coefficients {ck,s{f)}s£j<i(k) ™d reorder the indexes s G J'*(/c) as {sj}^Si^ so that 

|cA:,si(/)| > \Ck,S2{f)\ >■■■ \ck,sM)\ > ■■■ |Cfc,m(fc)(/)l, 
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and then define 

n{k) 

Gkif) ■■= ^{cfc,^^(/) - |cfc,s„(^j^j(/)|signcfc,s,(/)}Mfc,s.. 
i=i 

We prove that G S C{B'^g, Lq). For < r < oo, denote by V(A;)t- the quasi-normed space of all 
functions / E V(A;), equipped with the quasi-norrm L,-. Then by Lemma [3TT] qi. € C{B'^ Q,Y{k)p). 

Consider the sequence {ck,s{f)}s£jd(k) ^ element in i^^'^^ and let the operator : V(A;)p — > 

^"•-^^ be defined by 5 i-^ {o-sjsejd^k) if 5 € V{k)q and g = YlseJ''{k) o-sMk^s- Obviously, by 

Dk € C(S(A;)p, ^™''^^). For x = {a;fc,s}sgjd(fc) € we let the set {kj}^^ be ordered so 



that 

and define 



Pjl I > I > • • • I > • • • > Fj„(;.) I 

n(fc) 

^n(fc)(a;) := ^{xk, - k„(fc)+i|signxfc^)efc^.. 
i=i 

Temporarily denote by H the quasi- metric space of all x = {a^fc,s}sgjd(fc) € for which Xk = 

0,k ^ Q, for some subset Q C J'^ik) with |g| = n(A;). The quasi-metric of H is generated by the 
quasi-norm of £^^^\ By Lemma (|3.3p we have Pn{k) G C{£^,H). Consider the mapping i?M(fc) 
from H into S„(^)(M(A;)) defined by 

if X = {a;A;,s}sgjd(fc) G H and x^ = 0, k ^ Q, for some Q with jQI = n(/c). Since the family M(k) is 
bounded in Lq, it is easy to verify that RM{k) ^ C{H,Lq). We have 

Gk = RM{k) ° Pn{k) ° Dk°qk- 

Hence, € C{B^g, Lq) as the supercomposition of continuous operators. Since by Lemma [STT] 
Qkif) G C{B^^g,Lq), from it follows G € C{B^^„Lq). 

Notice that in the operator G, the quasi-interpolant Q^if) is the main non-adaptive linear 
part. Its adaptive non- linear part is a sum of continuous algorithms Gk for a continuous adaptive 
approximation of each component function qk{f) in the kth scale subspaces V(A;), k < k < k* . 

Let m be the number of the terms in the sum (I3.12p . Then, G{f) € Sm(M(/i;*)) and 

k* 

m = {2^+r-lf+ n{k). 

k=k+l 

Moreover, the number of sampled values defining G{f) does not exceed 

m' -.= {2^ + lf +{2^i + 2rY ^ n{k). 

k=k+l 



16 



Let us select k, k* and a sequence {n{k)}'^^j^_^-^. We define an integer k from the condition 

Ci2* <n< C22'^^, (3.14) 
where Ci,C2 are absolute constants which will be chosen below. 

Notice that under the hypotheses of Theorem ll.il we have < 6 < a. Further, we fix a number 
e satisfying the inequalities 

< e < {a-6)/6, (3.15) 



where S := d{l/p — l/q). An appropriate selection of k* and {n{k)}^_j^^^ 

k* ■= [e-^ log(An)] + a; + 1. 

and 



IS 



n 



(k) = [Xn2-<''~''^], k = k+l,k + 2,...,k* 



(3.16) 



(3.17) 



with a positive constant A. Here [a] denotes the integer part of the number o. It is easy to find 
constants Ci,C2 in (j3.14p and A in (j3.17p so that n{k) < m{k),k = k + l,...,k*, m < n and 
m' < n. Therefore, G is an n-sampling algorithm of the form ()1.4p with B = Sm(M(/c*)) and 
M(/c*) £ G 8lS a finite family. Let us give a upper bound for ||/ — S^{f)\\q. For a fixed number 
< r < mm{p, 1), we have by (12. Sh . 



(3.18) 



k* 

\\f-s^{f)\\i < \Mf)-GMfmi + Eii9^(/)iir 

k=k+l k>k* 

By ^M-^M and (IXTTD we have for all / G SB^g 

\\qk{f)\U « 2-("-^)^ A; = r + l,+2,... 



(3.19) 



Further, we will estimate \\qkif) — Gk{qkif))\\q for all / G SB^q and k = /c + l, k* . From Lemma 
I we get 

•m(k) 



1/9 



U=n(fc)+1 



< {n{k)}-'\\{ck,s{mp,k- 
By (IISD-dZj]), (l3l9D and (fCTj) we obtain for all / G S'S;^^ and A; = + 1, k* 



(3.20) 



kyk(/) - Gk{qk) 



m(k) 

E Cfc,.,(/)Mfe„ 

j=n(fc)+l 



m(k) 

E 

U=n(fc)+1 



1/9 



(3.21) 



« 2-'=/''{n(A:)}-iK,(/)}||p,fc « 2-"^2^'={n(A:)} 



-5 



17 



From ^J8\i by using IK2T^ . (f3lB . (f3ia - (i3Tril and the inequahty a > S, we derive that for all 
functions / G SB^g 

k* oo 

fc=fc+i fc=fc*+i 

k* oo 
^ ^-T"'52-T(a-<5)fc 2^'^('^^<^+'5s)(^^^) _)_ 2^'^('^~'^)^* 2^'^(°^'5)(*=^^*) 

k=k+l k=k*+l 

Summing up, we have proven that the constructed n-sampling algorithm G = S^{f) € C{B!^q, Lg) 
and is of the form (jl.4p with A = T,m(M.(k*)), and M(fc*) € ^ as a finite family for which the 
inequality (|5.5p holds true for the case p < q. U 



4 Lower bounds of Lq) 

To prove the lower bound Theorem ll.il we compare Vn{B^Q,Lq) with a related non-linear n-width 
which is defined on the basis of continuous algorithms in n-term approximation. 

Let X, Y be quasi-normed spaces and X is a linear subspace of Y . Let be a subset in X and 
$ = {vJfcjfcGA' ^ family of elements in Y . Denote by QiX) the set of all bounded families ^ dY 
whose intersection <1> nL with any finite dimensional subspace L in y is a finite set. We define the 
non-linear n-width {W,Y) by 

T^{W,Y) := inf inf sup ||/ - 5(/)||y. 

■fescy) SGC{x,y): 5(x)cE„(<i.) f^w 

Since all quasi- norms in a finite dimensional linear space are equivalent, we will drop "X" in the 
notation t^{W, Y) for the case where Y is finite dimensional. 

Denote by SX the unit ball in the quasi-normed space X. By definition we have 

MB^,e,Lg) > T^{SB^^s,Lg), (4.1) 

where we use the abbreviation: B := B^g. 

Lemma 4.1 Let the linear space L be equipped with two equivalent quasi-norms \\ ■ \\x o,nd \\ ■ \\y, 
W a subset of L. If {W,Y) > 0, we have 

T^+^(,W,X) < T^iW,Y)T^iSY,X). 

Proof. This lemma can be proven is a way similar to the proof of Lemma 4 in [15] . □ 
Lemma 4.2 Let < q < oo. Then we have for any positive integer n < m 

TniBZJ"^) > i(m-n-l)i/^. 
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Proof. We need the following inequality. If is a compact subset in the finite dimensional normed 
space Y, then we have [15] 

2Tn{W,Y) > bn{W,Y), (4.2) 

where the Bernstein n- width bn{W,Y) is defined by 

6„(W,y):= supsup{t>0: tSY n L^+i C W} 

Ln + l 

with the outer supremum taken over all (n + l)-dimensional linear manifolds Ln+i in Y. 
By definition we have 

Hence, by ()4.2p . Lemmas 13.31 and 14.11 we derive that 

1 = 6„_i(S-,C) < 2r^_i(i?-,C) 
< 2r„(i?-,^™)r„_„_i(i?-,C) < 2(m-n-l)-V«T:„(B-,^™). 

This proves the lemma. □ 

Theorem 4.1 Let < p,q, 9 < oo and a > 0. Then we have 



Proof. By (14. ip the theorem follows from the inequality 

T^{SB^^g,L,) » n^"/'^. (4.3) 

To prove (j4.3p we will need an additional inequality. Let Z is a subspace of the quasi-normed space 
Y and W a subset of the quasi-normed space X. If P : y ^ Z is a linear projection such that 
||^(/)||y < '^II/II'k('^ > 0) foi' every f gY, then it is easy to verify that 

T^iW,Y) > \~\^{W,Z). (4.4) 

Because of the inclusion U := SB^ g C SB^q, we have 

T^{SB^^e,L,) > T^{U,L,). (4.5) 

Fix an integer r with the condition a < min(2r, 2r — 1 + l/p, 2r). Let U{k) := {/ E V(A;) : ||/||oo < 
1}. For each / G V(A:), there holds the Bernstein inequality [llj . 

Il/lb^,, < C2"'=||/||oo, 

where C > does not depend on / and k. Hence, C^^2^°^^U (k) is a subset in U. This implies the 
inequality 

T^{U,L,) » 2-"V„^(^(A;),L,). (4.6) 
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Denote by V(/c)q the quasi-normed space of all functions / G V(A;), equipped with the quasi-norm 
Lg. Let Tfc be the bounded linear projector from Lg onto V(A;)q constructed in [11] such that 
< A'll/llq for every f £ Lg , where A' is an absolute constant. Therefore, by ()4.4p 

r^{Uik),Lg) » r^iU{k),V{k)g) = rn{U{k),V{k)g). (4.7) 

Observe that m := \J'^{k)\ = dimV(A:)g = (2'^ + 2r — 1)*^ x 2*. For a non-negative integer n, 
define m = m{n) from the condition 

n X 2* X m > 2n. (4.8) 

Consider the quasi-normed space of all sequences {as}seJ'*(fc)- Let the natural continuous linear 
one-to-one mapping 11 from \'{k)g onto £™ be defined by 

n(/) := {as}seJ''{k) 

if / G Y{k)g and / = E.ejd(fc) We have by ^-^1 ll/IU ^ \mf)h^ and ||/||, x 

2-'^'=/'?||n(/)||fm. Hence, we obtain by Lemma li^ 

T^iUik),Vik),) - 2-*/V„(BS,C) 

> 2-'^'=/'?(m-n- 1)1/'? » 1. 

Combining the last estimates and (|4.5|) ~ (|4.8|) completes the proof of (j4.3|) . □ 



5 Adaptive non-continuous sampling recovery 



In this section, we prove the asymptotic order of s„(-Bp g, M, Lg), r„(5i?pg)g and e„(S'i?pg)q in 
Theorem 11.11 

Let W and -B be subsets in Lg. For approximation of elements from W by B, the quantity 

£;(iy,i?), := sup mfJf-^Wg 

gives the worst case error of approximation. 

Let <1> = {'^k}k&K be a family of elements in Lg. The quantity of n-term approximation 
aniW, ^)g with regard to is defined by 

an{W,^)q := S(VF,S„(^'))g. 



Given a family B of subsets in Lg, we can consider the best approximation by B from B in 
terms of the quantity 

d{W,B)g:= mf^E{W,B)g. (5.1) 
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Notice the following useful identities 

aniW,^), = inf sup \\f-S{f)\\,. 

5:VF-i>S„($) f^W 

and 

d{W,B), = inf inf sup ||/ - 5^(/)||,. (5.2) 

The quantity d{W,B)g is called the entropy n-width (entropy number) e„(VF)q if B in (j5.ip is 
the family of all subsets B of Lg such that \B\ < 2". The non-linear n-width pn(W)q is defined 
only when Lg is a space of real- valued functions on a set 17, if in (j5.ip is the family of all subsets 
in Lg of pseudo-dimension at most n. 



From (|5.2p we have 



en{W)g = ^ inf inf sup ||/-5^(/)|| 



B\<2" SB-.W^B f^w "" -'-9' 



and 

Pn{W)g = inf inf sup |t/-S^(/)||,. 

dimp _B<n SB.W^B f(zw 

Therefore, we can take the last identities as alternative definitions of £n{W)q and pn{W)q. 

Theorem 5.1 Let p,q,9,a satisfy Condition (jl.Sp and a < 2r. Then for the d-variable Besov 
class SB^g, we can explicitly construct an n-sampling algorithm with B = S„(M) so that 

sup ||/-5f(/)||, X Sn{B^^g,M,Lg) X n-"/-^. (5.3) 



Proof. In \18\ Corollary2.3, Theorem 3.2] an n-sampling algorithm with B = S„(M) was 
explicitly constructed such that 

sup \\f-S^{f)\\g « n^-/'^'. 

feSB^,9 

This proves the upper bound of (15. 3p . 

The lower bound follows from the inequality s„(i?p M, Lg) > (7„(S'i?p g, M)g and the in- 
equality 

a„(Si?-„M), » n^"/'^. 
which was proven in jl8t Theorem 5.1] □ 

Theorem 5.2 Lei p,q,9,a satisfy Condition (jl.Sp . T/ien /or i/ie d-variable Besov class SB^g, 
there is the following asymptotic order 

rn{SB^^g)g X n-/'^. (5.4) 
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If in addition, a < 2r, we can explicitly construct a subset B in S„(M) having dimp(i?) < n, and 



a 



sampling recovery method of the form (jl.4p , such that 



sup ||/-5„^(/)||,«n-°/'^. (5.5) 



Proof. The inequality (j5.5p and therefore, the upper bound of (j5.4p was proven in |18[ Theorem 
3.1]. 

The lower bound follows from the inequality r„(S'i?pg)g > Pn{SB'^Q)q and the inequality 



p„(5R%), » n 



-a/d 



which was proven in [18t Theorem 5.3] □ 



Theorem 5.3 Let p,q,6,a satisfy Condition <\1.8\i . Then for the d-variable Besov class SB^g, 
there is the following asymptotic order 

en{SB-g), X n--/'^. (5.6) 

If in addition, a < 2r, we can explicitly construct a subset B in S„(M) having \B\ < 2'^ , and a 
sampling recovery method of the form (jl.4p , such that 

sup ||/-5„^(/)||,«n-°/'^. (5.7) 

Proof. The inequality (|5.7p and therefore, the upper bound of (|5.6p was proven in [181 Theorem 
4.1]. 

The lower bound follows from the inequality e„(5i3p g)^ > e„(5i?p g)g and the inequality 

-a/d 



en{SBpQ)g » 



n 



which was proven in [18^ Theorem 5.5] □ 
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